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AN AUTOMATED PROCEDURE FOR CALCULATING SYSTEM MATRICES FROM 
PERTURBATION DATA GENERATED BY AN EAI PACER 100 HYBRID 


COMPUTER SYSTEM 

by Edward J. Milner and Susan M, Krosel 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

Most control system design techniques assume having available a 
linear representation of the system to be controlled. Although obtain- 
ing a linearized model of a complex, non-linear system is usually tedi- 
ous and time consuming, this task can be made less so by using the 
automated procedure described in this report. 

Techniques are presented for determining the elements of the A, 

B, C, and D state variable matrices for systems simulated on an EAI 
Pacer 100 hybrid computer system. The techniques are generalized to 
allow use with many different simulation models. An automated proce- 
dure systematically generates disturbance data necessary to linearize 
the simulation model and stores these data on an auxiliary TEKTRONIX 
model 4922 floppy disk. 

A separate digital program verifie he floppy disk data, allows 
for data modification, if necessary, cai., ■ -,.v cs approximations to the 
partial derivatives forming the elements of the system matrices, and 
prints the matrices appropriately labeled. 

Use of the floppy disk as a storage medium permits the data gen- 
eration process on the hybrid computer to be separate from the matrix 
element calculation process. 


INTRODUCTION 

The intricate nature of many dynamic systems necessarily requires 
rather complex simulation models to represent performance over their 
entire operating range. One is likely to find the model to include a 
combination of linear and nonlinear algebraic and differential equa- 
tions and maps. Although such a model is useful for predicting system 
performance over a wide range of operating conditions, generally, its 
nonlinear characteristics hinder it from being well-suited for control 
analyses. A basic assumption of most control design techniques is hav- 
ing available a linear representation of the system to be controlled. 



By considering operation in tlie vicinity of a particular oper~ 
atlng point I state variable techniques can be employed to obtain a 
linear approximation of the nonlinear system. Doing so, the system 
equations can be written in vector notation as: 


X “ Ax + Bu 

(1) 

y « Cx + Du 

(2) 

where x Is a vector of the system time-varying states, x is a 
vector of the time rate of change of the system states, u is a vec- 
tor of system inputs, and y is a vector of system outputs. Matri- 
ces A, B, C, and D are characterized by: 

n - 

aij - ^ 

C3) 


(4) 


(5) 

d 

“ij 3uj 

(6) 

One linearization technique is to obtain finite difference approxi- 
mations for the partial derivatives which comprise the elements of 
the system matrices. To generate the partial derivatives, the steady' 


state level of each state variable and system input is Independently 
stepped In both a positive and negative direction while holding the 
other state variables and system Inputs fixed. Corresponding values 
of the state variable derivatives and system outputs are recorded for 
each step change. By calculating the deviation of each state vari- 
able derivative and output from its steady-state value and by form- 
ing appropriate quotients the matrix elements are formed. 

This report describes an automated technique for determining 
Uie elements of the A, B, C, and D state variable matrices for sys- 
tems that are simulated on an EAI Pacer 100 hybrid computer system. 
The technique is generalized to allow its use with many different 
simulation models. The data taking process is initiated by manually 
setting a sense switch at the Pacer 100 control console. The order 
in which the variablee are to be perturbed along with other appropri- 
ate data are input from cards. Parameter values are transmitted from 
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the analog computer to the digital computer by means of analog-to- 
digital converters (ADC's). Perturbation data generated on the Pacer 
100 system are recorded and stored on a TEKTRONIX model 4922 floppy 
disk. The data are later transferred from the floppy disk to an IPM 
TSS/360 computer where the matrix calculations are performed. The 
matrix calculation is contained in the program MATCALC which can 
generate matrices for systems having an arbitrary maximum of 20 
states, 10 inputs, and 20 outputs. 

The following sections provide details of the apparatus, the 
digital program, the perturbation process, the floppy disk data for- 
mat, and the matrix calculation procedure. An example is provided to 
Illustrate the procedure. 


APPARATUS 

The data taking technique, described in this section, applies to 
nonlinear system simulations impleuented on an EAI Pacer 100 hybrid 
computer system. The technique could, of course, be modified for 
other hybrid computer systems. At vhe NASA Lewis Research Center, 
the Pacer 100 system includes the Pacnr digital computer with 32K of 
core, a TEKTRONIX model 40.10 teri>:lnal with auxiliary TEKTRONIX model 
4922 floppy disk unii , two model 6bl analog computers, and an inter- 
face to allow communication between the computers. That Interface 
consists of 48 analog-to-digltal converters (ADC's) and 24 dlgital- 
to-analog converters (DAC's). 

Analog Setup Requirements 

It is assumed that the integrations associated with the system 
dynamics are performed by the Integrators on one or both of the ana- 
log computers. The integrator outputs then correspond to the system 
state variables. The perturbation program requires that each state 
variable, state derivative, input variable, and output variable be 
transmitted from the analog computer to the digital computer through 
the interface. Hence, each parameter of interest must be connected 
to an ADC. For the purposes of this report, parameters formulated 
on the "primary" analog console are transmitted to the digital com- 
puter via a block of ADC's beginning with ADC 1 and parameters form- 
ulated on the "secondary" analog console are transmitted via a block 
of ADC's beginning with ADC 33. 

Digital Program Description 

Appendix A contains FORTRAN listings of the major parts of a 
typical simulation program which incorporates the necessary pertur- 
bation logic. 

The program MAIN is designed to control the execution of the 
major parts of the program. As the execution of each part is 
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completed, control Is returned to MAIN. The analog computer eetup 
and checkout takes place during the execution of subroutine SET. 

This includes setting and checking analog pots, digitally controlled 
attenuators (DCA's), and digitally controlled function generators 
(DCFG's). Subroutine STATIC provides for the reading of all ampli- 
fier outputs at any desired operating point. Subroutine LOOP is the 
heart of the simulation insofar as this is the area where the dynam- 
ics of the system being simulated ere updated. By continually 
cycling through component maps and table lookups, current values of 
the system state derivatives are obtained. In general, LOOP will re- 
ceive information from the analog computer and output information to 
the analog computer through the interface. 

The perturbation program consists of three primary digital sub- 
routines - DATOUT, SGNFIX, and AVG. The function served by sub- 
routine DATOUT is to set Che Integrator initial condition DCA's to 
correspond to the steady-state operating point being considered and 
to systematically perturb the state variables and the system Inputs 
one at a time in both a positive and negative sense. DCA's are used 
because they can be set Ly the digital computer without changing 
modes on the analog computer. Subroutine SGNFIX is used to guaran- 
tee chat the positive algebraic sign of each parameter is recorded 
on floppy disk. For example, consider a variable, x, whose range 
spans the set of all real numbers. Hence, x may have both positive 
and negative real values. However, the analog computer representa- 
tion can be either the variable (+x) or the variable (-x) depending 
on Che sign of the reference voltage used. If (-x) were the analog 
computer variable, subroutine SGNFIX converts it to (■♦■x). For each 
perturbation, subroutine AVG scans the parameters of interest form- 
ing 100 point averages for each of them. These data are then re- 
corded on floppy disk. 

Formatting the output data and recording it on floppy disk li 
accomplished in subroutines ;.nDL.i -aiu* DIS!’ » respectively. Tl.c 9 <‘e 
purpose of subroutine is to control the recording of the data 

on floppy disk. To record on floppy disk, the floppy disk unit is 
plugged into the disk interface card on the TEKTRONIX 4010 CRT term- 
inal. The write statements in the program are the same as if the 
data were being written on the CRT terminal. Data is then automa- 
tically written on the CRT terminal and recorded on the floppy disk. 
Note that everything written on the CRT terminal is also recorded 
on the floppy disk. 


Illustrative Example 

Consider the following system of four uncoupled dynamic systems 
linked through output relationships: 
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ST81 - J(INPUTl-ST81)dt 

where 

1 INPUTl 1 

i_2 

and 

|ST8l| <2 

(7) 

ST82 -J(INPUT2-ST82)dt 

where 

1 INPLT2 1 

12. 

5 and 

|ST82| £2.5 

(8) 

ST83 -/(INPUT3-ST83)dt 

where 

|INPUT3| 

^5 

and 

|ST83| < 5 

(9) 

ST84 -/(lNPUT4-ST84)dt 

where 

|INPUT4| 

1 ^ 

and 

|ST84| <_4 

(10) 


An analog diagram representing a typical system such as these Is 
presented In Figure 1. Equations (7) and (8) are to be programmed on 
analog console 1. DCA 45 Is to be used as the initial condition pot 
associated with ST81 and DCA 41 is to be ured to set the perturbation 
amplitude of INPUTl. Likewise, DCA's 75 and 71 shall be used with 
ST82 and 1NPUT2, respectively. 

Correspondingly, Eqs. (9) and (10) are to be programmed on ana- 
log console 2. However, on this console DCA's 75, 71, 45, and 41 
shall be used with ST83, INPUT3, ST84, and INPUT4, respectively. 

In all cases a 2 percent perturbation amplitude is deemed suffi- 
cient. 

The solution of this example problem will be followed through in 
the discussion which follows. 


PROCEDURE 

The data taking process is activated by setting sense switch 4 
and then exiting from the main system calculation subroutine (called 
LOOP) by setting sense switch 1. Control is then transferred to 
DATOUT. Hereafter, DATOUT controls reentry into u(X)P. An internal 
check is made that the analog computer was in the operate mode at the 
time the data taking process was activated. Referring to Figure 2, 
a card is then read in (313) format giving the number of state vari- 
ables, the number of system inputs, and the number of system outputs 
formulated on the primary analog console only. A series of cards 
punched in (F5.1, 13, S7) format are then read. Each card contains 
information characterizing one of the parameters of interest. For 
each state variable, the algebraic sign of the integrator output, 
the Initial condition pot number, and the proportion of perturbation 
amplitude desired are required. For each system input, its algebraic 
sign. Its perturbation pot number, and the proportion of perturbation 
amplitude desired are required. For each system output only its 
algebraic sign (as previously discussed in the explanation of sub- 
routine SGNFIX) is required. 
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The next card read la a card In (313) format giving the number of 
state variables, the number of system inputs, and the number of system 
outputs formulated on the secondary analog console. Simllarlly, 
another series of cards In (F3.1, 13, S7) format are then read charac- 
terizing the parameters of Interest on the secondary analog console as 
was just done for the primary console. Table 1 contains a listing of 
the cards used for the example and Figure 2 shows the data deck setup. 

The order In which these data cards are read Is very Important. 

The parameter order determined by the cards must correspond with the 
parameter order determined by the analog ADC's used. State variable 
data Is followed by system input data, which, in turn, Is followed by 
system output data. Also, the order In which the parameters appear 
within each of the first two groups is important Insofar as the order 
in which the cards are read determines the order In which the pertur- 
bations will occur. Note that the order of parameter perturbations 
is in Itself arbitrary. Limitations on this order are imposed only 
by the analog console on which the integrator for a state variable or 
system input is located. As this procedure is defined, all state var- 
iables and system inputs on the primary analog console will be per- 
turbed first. Then those on the second analog console will be per- 
turbed. Table 11 gives the order in which the state variables and 
system inputs of the example are perturbed. Also, given is the ana- 
log console for each parameter and its associated ADC channel. 

Following the reading of these data cards, the digital program 
proceeds to monitor the steady-state integrator outputs and sets DCA's 
corresponding to the integrator initial condition values. After the 
steady-state operating point is established, the analog computer is 
automatically put in 1C mode where it remains during the perturbation 
and data taking process. The data taking process begins with a base- 
line reading of the staady-state values of the parameters of interest 
and the recording of them on floppy disk. The set of baseline data 
is designated as reading 1. Following the taking of this baseline 
data, each state variable and system input Is stepped sequentially in 
both a positive and negative sense. Table 111 identifies the pertur- 
bations associated with each reading for the example. 

For each perturbation, the data written on floppy disk are aver- 
ages of the parameters of Interest. The ADC's are read s«. quent ially 
100 times with the corresponding values summed and averaged. Since a 
parameter may be represented on the analog computer with either posi- 
tive or negative algebraic sign, each average is multiplied by either 
(+1) or (-1), depending on the parameter sign. Thus, the data recorded 
on floppy disk corresponds to the positive algebraic sign of each 
parameter. 

Data associated with each perturbation is recorded on floppy disk 
and each set of data on floppy disk is assigned its own identifying 
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reading number. Reading numbera are Incremented as each set of data 
Is recorded on floppy disk. Data is recorded on floppy disk according 
to the following format. The reading number is recorded twice In 
(215) format. The reason for this redundancy Is so that it may be 
used as a check when reading the floppy disk. The parameter averages 
are recorded on floppy disk as integers (with 3276?2 ^q representing 
full scale on the 32K Pacer) at the rate of 18 values per line to- 
gether with a line checksum according to the format (1816,17). The 
checksum value is the sum of all the integers output on that particu- 
lar line. If the number of averages exceeds 18 but Is not evenly di- 
visible by 18, the second last line of output will be the n remain- 
ing averages recorded in (nI6) format and the final line will be the 
corresponding checksum In (17) format. If the number of averages Is 
less than 18, the n averages will be recorded In (nit) format with 
a second line containing the checksum in (17) format. An example of 
the floppy disk data format la given in Table IV for the test case. 

The final characters required on each line are a "line feed" and 
"X-off". This ending sequence is necessary to insure proper data 
transfer to the IBM TSS/JbO computer. Because this is not the normal 
line termination, the contents of the digital computer core location 
controlling line termluatlcn must be changed. The address of this 
core location Is obtained by taking the address of .FOR from the mem- 
ory map (for the listing given li. Appendix A, this address Is 17234g) 
and adding 213g. The address displacement value, 213g, depends on 
the version of the .FOR routine used. Version i:o , December 1975 Is 
in use on the Pacer 100 system at Lewis. After loading the program, 
this memory cell Is opened and Its contents set to zero by the user. 

The entire procedure is repeated for each perturbation. When 
the data taking process Is completed, an appropriate message Is out- 
put on the line printer. At this point the system may be returned to 
normal operation by resetting sense switches 1 and 4. 

MATRIX CALCULATION TECHNIQUE 

This section describes a matrix calculation technique which uses 
the IBM TSS/360 computer with data Input provided by the floppy disk. 
The matrix calculation is contained in the program MATCALC. This Is 
a generalized version of a program developed at Lewis for matrix gen- 
eration from data obtained from a hybrid simulation of the FlOO engine. 

The MAT(]ALC program has three basic sections: data input, data 

verification, and matrix generation The data Input section uses a 
software checksum to insure that tnc floppy disk data was correctly 
transmitted to the IBM TSS/360 computer. After the data has been cor- 
rectly transmitted, the data Input section unscales It and stores this 
unsealed data in an Internal program matrix for use in the other two 
program sections. The data verification section determines: (1) If 

only one state or Input was perturbed for each reading and (2) if the 
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parturbaclon size of each state or Input was at least IX of the state's 
or Input's steady~stat« value. The matrix generation section calculates 
the approximations to the partial derivatives forming the elements of 
the A, B, C. and D matrices and prints out the matrices In their final 
form. A listing oI the MATCALC progr;<« Is given In Appendix B. 

Transfer of data from the floppy disk to the 360 computer Is con- 
trolled by a proedef called READ IN. A listing of this proedef Is In- 
cluded In Appendix B and the procedure for use is given in Appendix D. 
The proedef, READIN creates the vs-dataset, 'DTRK' which contains the 
floppy disk data. The data on disk Is a series of readings where each 
reading Is a list of values of the states, state derivatives, inputs, 
and outputs. The required format of each reading was discussed in the 
preceding section. There should be 2n’^l readings on disk for each 
set of matrices, where n Is equal to the number of states plus the 
number of Inputs. The first reading Is the 'base run' which is a read- 
ing of the unperturbed, steady-state values of the states, etc. (NOTE: 
The MATQALC program terminates if there Is no 'base run'.) Correspond- 
ing to each reading is a disturbance number. The 'base run' must be 
assigned the disturbance number of 1. Since each state and input Is 
perturbed in both a plus and minus direction, there are two readings 
associated with each state and input. These readings should be as- 
signed disturbance numbers of m and ml where m is an even in- 
teger in the range 2 to 2n. 

The readings contained in dataset, 'DTRK', are used by the data 
input section of MATCAlC to form an internal program matrix, DATA. The 
rows of the DATA matrix are the readings in the same order as they ap- 
peared in the floppy disk data and in dataset 'DTRK'. Each row con- 
tains .he unsealed values of the corresponding reading (i.e., the nth 
row contains the unsealed values in the nth reading). The order of 
the values is the sane for the reading and the row. Table V lists the 
DATA matrix for the c^xample. 

The rows of the DATA matrix are used In the data verification and 
the matrix generation sections of MATCALC. The INDX array, also formed 
in the data input section, is used in referencing individual rows of 
the DATA matrix. Table V also Includes the INDX array for the example. 
This referencing Is accomplished by: 

Row No. of DATA matrix - INDX (DISTURBANCE NO.) 

If the Row No. is zero, the reading data associated with that disturb- 
ance number was not contained in the floppy disk data (i.e., the data 
for a plus or minus perturbation of a particular state or Input was not 
recorded) . 

The NAME array associates each state or Input with two disturbance 
numbers: ore for the plus perturbation and the other for the minus 
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perturbation. Thla array contains the namea of the state varlablea 
and the ayatem inputa In the order In vrhlch they are perturbed (see 
Table ll)* The data for the NAME array la contained In the Input 
dataset* 'QINFO*. 'QINFO' la a va-dataset created by the user and 
contains data describing the contents of the floppy disk and speci- 
fying the structure of the A* B* C* and D matrices. The format of 
'QINFO' Is given In Appendix C and Table VI lists X'\e 'QINFO' dataset 
for thn example. 

Appendix D contains a detailed procedure for procesaing of per- 
turbation data which has been stored on a floppy disk. Instructions 
for operating the TEKTRONIX 4010 CRT terminal and TEKTRONIX 4922 
floppy disk drive unit are Included. IBM TSS/360 commands for reading 
the floppy disk data Into the IBM TSS/360 computer and for performing 
the matrix calculations are provided. Also In Appendix D Is the pro- 
cedure to follow should checksum errors occur. 

Table VII shows the resulting matrix listing for the example. As 
the results show, the slmuK<red system consisted of four uncoupled dy- 
namic systems linked only ^bvough the output relationships. The lin- 


earized system equations are: 

x^ - -0.9938 x^ + 0.9922 u^ (11) 
X2 - -0.9918 X2 0.9888 U2 (12) 
X3 - -0.9939 X3 + 0.9878 U3 (13) 
X4 - -0.9923 X4 + 0.9861 U4 (14) 
yj - 1.003 xj^ - 0.001022 X3 (15) 
y2 * -0.004666 X3 1.002 X2 + 0.001022 X3 + 0.001019 U3 (16) 
yj - 0.999 x^ -0.001546 x^ (17) 
y^ - 1.002 x^ 


( 18 ) 
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Program Conaideraclona 

Some modification of the perturbation program would be neceaaary 
if one of the atate varlablea or Inputs were at a limiting value 
while perturbation data were being taken. Because of that limit, 
that parameter could not be perturbed in both a positive and negative 
sense. Hence, proper logic would have to be included in subroutine 
DATOUT to restrict the perturbation of that parameter to the direction 
which would not violate the limiting value. Likewise, if the problem 
sixe were such that only one analog console were required, subroutine 
DATOUT could be easily modified to eliminate parallel operationa ad- 
dressed to the secondary analog console. 

Use of a aeparate digital computer (the IBM TSS/360 computer) al- 
lows for the processing of the perturbation data and analysia of the 
reaulting matrices to be done without interfering with use of the 
hybrid computer system. Calculation of the system matrices could have 
been done on the Pacer digital computer with the perturbation data 
stored on the Pacer's associated disk. 

The format specified for data output to disk is arbitrary and can 
be changed if desired. The DISK subroutine of the perturbation program 
and the data input section of the matrix calculation program should be 
modified for the desired format. A new checksum error check should be 
developed to insure correct data transfer from disk. 

CONCLUDING REMARKS 

The A, B, C, and D state variable matrices used in controls design 
can be readily determined using the procedure described in this report. 
Data required to obtain a linear representation of the system to be 
controlled are systematically gathered using the techniques presented 
in the text. Disturbances are introduced in the nonlinear system one 
at a time. The user may determine the level of each disturbance and 
the sequential order in which each is to be introduced. For each dis- 
turbance, corresponding levels of essential ystem parameters are auto- 
matically recorded on floppy disk. 

Use of the floppy disk %s a storage medium permits the data gener- 
ation process on the hybrid computer to be separate from the matrix cal- 
culation process. This floppy disk data is input to a separate digital 
program which contains the matrix calculation process. A software 
checksum is used to insure correct data transmission from the floppy 
disk. This program also checks the data for proper disturbance levels 
and correct number of disturbances. Appropriate messages are output 
if errors are detected, thus allowing data modification, if desired. 

The partial derivatives forming the elements of the system matrices are 
approximated by finite difference calculations. The system matrices, 
appropriately labeled, are printed out by this program. 
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Ub« of this procodura on a dynamic ayatem almulatod on a hybrid 
computar ganarataa automatically the data raqulred for linearization 
of the ayatem and detarmlnea the matrlcaa which repraaant the ayatcm. 
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ORUILNAL PAliK lb 

OF POOH Ql’AU’l*V 


NATCAU'* rUKIAH AMO 'UAOIM' rAOCDIP 

bfH^PAttnk AbUrtAX 

{.tAihSIOV AAAXlJOt, l«0l|AI) 

DlRKNSlbR IVAL(IP). ti»rA(«1,T0| 

OllfcHlIUk AbAltdOl. »)bSV(?0|, NDMVdO) .kOlClOl 
oiftmiON ’iM/O'.ATisrtioi 

blltraiXOIl A (iO,20),b(/}, 10) ,C(20,20) ,0(X0. 10) 
DlAl.liXON MC>'](tO),Aba(JO) .H««7(20) ,kAA)(20) 

DAtA bkTA/A2f')*l).0/ 

OAIA .LAXT/0/, tP0I/H«0/ 
bATA PTaST/.'M/, AtPO»/0/ 

XHTfll* VA. ,Nb-i,AbO 

Dal A UbV*0*«H /, ADO/IOaAH / 

>KAb(H,3001) NS. N;, NT 
4001 PCaPA! (JTi) 

NnOAOI • 2AnS*NC*«T 
NA'jkS < 2* (MS'NC)* 1 
NI«T • kS»«C 
N'.xze • as 
aXNi'ii; • Nc 

k '.A b( 1,0002) (MAie(I) .I«1,btOT) 

0002 PokPAT (lOA.) 

kSAOd.kOOl) (SP(l) ;N0S) 

AOO} PulNAr ((< (10P4. 2/} ) 

kXAriA.OOO/) |IIAN2(I) ,I*1,NS) 
kXAb(1.*002) INAnl (i) .!• 1, NT) 


PLukk) DXSa AlAb IN AND 'IkSCALINQ OP DATA 
••ktAD XN STAPT 

x-c 

1 I«JO 

k2Abl%,k00t,tk0-2,tkP*2| kA, as 
0001 POkNAT <2XS) 

• •DlATOIlbaNCe NONbkh t.HA:k 
IP <NA ,ku. kk) 00 TO I 
NAITa (6.4100) NA.Nk 

6100 POkNAT r»*,*NA VALUE {'.n,*) IS NOT EOOAt TO NB VAtOl (', X 1, •)’./. * 
A61,'ENTIl.X COkkECT JISTUABANCl VALUE (IS) POkNAT') 
ktAMT.N 101) NA 
7101 P06NAT (IS) 










n n 




2) 


rtii TO J 

i KMTt (t, 700^1 

7002 »jrf»T or mm <» 't kjoal to rirrrTtr no. or ttisTurBAi^ts',/. - 

*Sl,*UiT»r 1 TO STOr, Urt'ILT TO kONTtHOt*/) 

K.t (7.700J) K^'isr 
700J rOM«T (II) 

ir (KTI : T .tv. •) STor 
r.o TO SoO 

I ir (i «DI (««) . 10 . 0) 10 TO « 

•*:tt (' , 100 1) Mk . IHI.I («*) 

1001 ro»r*T ibx.'iiiar (', u,») pruvrouni srr tj seo»,i»./.* 

• sx,*i>rr*ULT TO mri/ict, i to stnr«/) 

HAD (7,7001) KTESr 
If (NissT .ro. 1) tiff 

c 

C ••mi,* - CHiTTHTr AHr SIQUXtIC* A'JWBtrS 
It IhDMNA) • I 
ICBT • 1 
J • 1 
11‘ia • iH 
5 JA • J* 17 

ir (LCNT .tu. 1 .AMO. HEND .OT. NWOlOS) TO T3 10 
ir (J . >T. Ai>7ro&) GO TO 10 
rtAa (0,%002) (IVAk ( 7 ) , K>1. HlltD) , icr 
?002 FJ'TAI )1(»i6,l7) 

ICil • U«T»1 
GO TO 17 
10 .lA ■ Nixoros 
peso • JA-J»i 

rf.AS (‘>,S0J2) (IVAl, (A) , r> 1, MrNb) 
l>t*E (O.bOOi) KK 
5001 70M>AT (17) 

LONT • LCir*; 

••cHXCRS'll CIIECA 
15 ISU7 ■ 0 

liO 15 Kd.NLNb 
U 15UPI • tSII1*IVAL(K) 

ir (I7U1 .Nl. ICK) GO TO 50 

100 HAST • icr 

hUTTE (7,5001) SA,1A 
If (MFIIO .NE. IP) 30 TO 17 
A*tTE (7,5002) (IVAL(A) ,r«1,*' n),ICA 

GO ro 20 

17 hRITE (7,5002) ( I VA L (F) , K ■ 1 , MEND) 

7PITE (7,500)) ICR 
C 

(. ••UMS^ALIRG OP SA5PU0 VAPIANLES 
20 P* 0 

bo 25 M»J.JA 
R » R • 1 
ri • IVAl(K) 


- i 


l 


i 


n#>noonoonn 


.1 



2 * 


YS » Tl/J2‘>6*. 
a b»rm,») • fs'srjii) 

J ■ JAM 

ir (j .u. NrioiDS) ao ro % 

ir (1 .LT. NtailSI YO TY 1 
MMTC 

»001 (5l,MAY»II DATA PCAO IN AND UNSCALIHY CON' CT’'/,- 

A'^K.'OATASAT UN UNIT 9 CONTAINS (CONNACTtO IP NIC) PLOPPT OISN DATA*) 

80 TO SOO 
C 

C ••CrtLCRSUH PNNON NOUTINE 

)0 NAITC (N.6I0I) NA. I, L:NT 

6101 ACkSAT (• CHfcCKSOP AM0N«,1I,' (UN '.II,' S19 '.I>« ' LIHI! '.INI 
U>m (6,6102) HAST, ICR, ISilN 

6102 rOPNAT (• CHECP309S l PREVIOUS ■ *,T7, • .YEAD • ‘.IT,* CALC • *,I7) 

NAITS (6,6101) (IVAL(K) , K« 1 , NANO) 

6101 POARAT ( ' THE VALUES APE • , /, I I , 1SI6 ,/) 

200 NAITZ (0,6104) 

6104 POPRaT (• ENTAP CJRAECTEO LINE*) 

IP (NANO .20. IS) 30 TO 6010 
PAAD (7, 6002) (IVAL (K) , «■ 1 , NEMO) 

AAAO (7.LO0J) ICP 

80 TO 6920 

6010 |it.AD (7,5002) ( I VA L (P) , K « 1 , NEND) , ICP 
6020 isun • U 

I'j 55 K>1,NA40 
55 ISUR ■ XSUR*IVAL(P) 

IF (tSUR .FfJ, ICP) 00 T9 lOU 
4MTE (6,6105) I3UR ,ICP 

6105 P8PRAT (• ChECPSURS DIPPER, CALC • *,I7,* READ • *,17) 

HKIIE (6,6106) 

6)06 PCJRAI (• ENTER 1 TO PE-ENTEP , 2 TO 3T0P,0tPAULT TO CONTINUE*) 

PAAD (7,700J) HTEST 
NTESI » NTEST*1 
«0 TO ( 100,200, 300) ,NIEST 
300 STOP 
C 

C 

500 CONTINUE 

ir DESIRE A LI3TIN3 OP DATA RATPIE NHICH RON CONTAINS ONSCALBO OISTORBANCl DATA, 
INSERT RRITE STATERANTS HERE 


• • **•••• 


INITIAUIATION POP DISf'IRLANCE VERIPICATION 

.4NPUN « IHOI (1) 

IP (NBRUN .NE. 0) 00 Ti 544 
4R1TE (6,6001) 


27 


6002 roanAT (fti.'No i>asK Dara.TEtN ihatc* ) 

JTOP 

6I*D|H,400%) (N0ILT(K| .a«1,aTOT) 
ar«b((t,ooo%) (HOatv (at ,a<i, Mf) 
rekoiB.oooS) (MOBsv(ai 
FMb (a,800S) |Mcoua).K>i,aruTi 
ttoos rortAT(ioii) 
c 

1,0 510 a»1,NT0T 
K3 ■ !IDELT(M) 

510 ATfcST(N) ■ Sr(ND)/1000. 

c ..••••••• 

c 

c LisToaitANce check and treairiCArioN 

c 

DO 700 M«1,NT0T 
l•k•M•NSI^t 
Ii33 * NANE(N) 

DO 600 1*1,2 
NC*0 

NSX • 2»A*1-1 
M ■ IN07 (NM) 

If (NI .EO. 0) NX > NBNJN 
to nOO X*1,M0T 
.•J « NOaLT(X) 

TP (ASS (D/ TA (NBBON, 1D|) .LT. ATESTUM SO TO 199 
TEST • (D «TA (NX, NQ1 'DAT A (MBBUM, NDI) /data (NB»UN,NO| 

IP (AH'.riST) .LT. PTESTI SO TO NOO 
112 XP (I .Nb. Ml 00 TO 331 
NC • X 
r,0 TO 400 
313 Nil ■ NAIEdl 
IKloK • NrfcO(t*1 

NPITE (6,201) NI,HflI, Nil , hi 3, DATA (N1,HD) , DATA (NBBOM,ND) 

201 fORIAT (1H ,'StwUENCL ',I2,> INDI •,I2.2N. ,A4,' IS FENTOPBEO IN ADDITION TO «,A4,/, 
AN|,'VALUI •',013. 5,* BASE VALUE • ',011.5) 

*aXTE (e,201| NI 

203 POR'IAT (* ENTER 1 TO USE BASE INSTEAD OP SEO ',12,' OTHEPNISE DEPAULT'I 
hEAl (7,7001) NTESr 
IP (NTES1 .NE. 1) 00 TO 400 
XNDA(NI>I| ■ 0 
NPITE(6,204| NI 

204 POPHAT (51, 'BASE IS bEINS USED POR SEQ ',12,/) 

SO TO 400 

399 TEST * OATA (NX ,ND) -DAT A ( NBRUN,ND) 

XP (ABS(TEST) .ST. 2»ATEST(I)) SO TO 332 

400 CONTINUE 
XP (NC . Eg. N) SO TO 600 
If (INOX(NRX) .Eg. 0) SO TO 590 
NEROP » NERni'*1 
N D • N D E LT ( i> ) 

*RITE ({,202) NI,NRI,N3 J.DATA (MI,N.,) ,DATA |NBR'IN,ND| 
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202 rOf«»T (Irt ,'SkguetlCl •,I2,», moi M2,*, OISTOIIBtHCI to* OB 

• •,011.5, • BkSC VILUC • ',011.51 

r.O TO 600 

590 MBITS (6,209) NRI,l«Sl 

209 rOBlUr (• IBDX •,I2,*, BO DI5TUIB1MCE BOB »,»4,/,- 
USIBI B5Sr BUN*) 

600 CONTINUE 
700 CONTINUE 


PK7*U CNLCULtTlCN 
J»« 1 

10 600 N«1,HT0T 
NCL ■ NCOL(N) 

NN ■ N'NSIZE 
NS) « NNHE(NrL) 

.XMORO • NDELT(NCL) 

NTbO • 2*NCL 
NIOHl ■ INDX(NTMO) 

NB0M2 ■ INDX(NTH0*1) 
ir (NB0B1 .Eg. 0)N6ON1 ■ NBRUN 
ir (NR0N2 .EO. 0) NB0H2 ■ H6BUN 
607 PNiJNN ■ DATA (NBOHl,NWOBr<) -DATA (NB0M2,NH0BD) 
ir (ABS(DNOHk) .or. .001) 00 TO 609 
IB (NB0M1 .E3. NbRUN) 00 TO 60S 
Nt ■ HI0N1 

NRITE (6,205) NSJ,NI 

205 MR1AT (1H ,AU,« DISTURBANCE lObNTICAI, SETTIN3 SEOUENCE ',I2,' TO BASE') 
N10N1 ■ NBRUN 
GO TO 607 

60b NRITE (6,206) NS3 

206 hORNATC NO ',AU,' DISTURBANCES, NATRIX COLUMN ZEROED') 

LSOnk » 1.0E50 

C •••CALCULATION OB A e B NATRIX COLUMNS 

609 LO 620 K<1,NSIZE 
R3Un ■ NDRIV(K) 

IB (K .31. NSIZE) 30 TO 610 

A(K,K) ■ (DATMNB3H 1,NDUN) >DArA (UR0N2,NDUN) ) /O'lONN 
GO TO 620 

610 B(K,NN) < (DATA (NR3W1, NOUN) -DATA (NR0H2,NDUN)) /DNONN 
620 CONTINUE 

C •••CALCUALATION OB C 6 D NATRIX COLUHNS 
CO 610 K«1, NSIZE 
KJ'JN • NORSV(K) 

IB (N .r.T. NSIZE) 30 TO 625 

C(K,N) « (DATA (NP3U 1, NOUN) -DATA (Nt>0k2,NDUN) ) /DNONN 
GO TO 6)0 

625 C(K,KN) « (DATA (NR0H1,NDUN) -DATA (NR0M2, NOUN)) /DNONN 


I 


/ 


o o o n n 
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t)0 C >NTttiUE 

J»J» J 

ik» > NCOUlil 
:*(ii .f!T. aj TO 710 

l('(» .OT. W3) GO ro 70S 
hJ'i(J) « N4ll(VkP) 

GO TO BOO 
705 J» 3 

710 • tiSPMVtR) 

BOO CONTIKUE 




I MkTUUT ur IATFIClS 
10 COLS n*x PER Like 
SC - Ri^ipiir 
«3ITi (b,k*00) 

6000 rO.-SAT (bX.'TUE A kATRUV) 

IK (SS .Lt. 10) GO TO o950 

MCOLI'R'l/^ 

li<:OL2>NCOL1« 1 

Aim (0,69^0) (HDG(L) ,L>1, JO) 

6960 ropfiArr •,9 x,joa4) 

LI 6906 t«1,KS 

6906 H,T1TF (6,6J01) h*R2 ( I ) , ( A ( t , J) , J> 1 , KCOL1 ) 
kllTE (6,» 960) (HDG (L) , L> 11 ,60) 

DO 6907 1-1, RS 

6907 90ITE (6,6901) kAR2 (I) , ( A (I , J) , .)■ NCOL2, NS) 
GO TO 6955 

6950 WH'TE (6,6960) (HDG (L) , L< 1 , 10) 

CO 6908 1-1, NS 

6906 kniTL (6,6901) kAN2 (I ) , ( A ( ( , J) , J> 1, KS) 

6901 fORNAT (• • ,A4,5E, 10G12.9/) 

6955 HHITL (6,6902) 

6902 rOKNAT (51, 'THE B NATBIX'/) 

Wt'tTE (6,6960) (AOG(L) ,L-1, JO) 

C'J 6911 1-1, NS 

6911 NNITE (6,6901) NAN2 ( I) , (0 (I , J) , J- 1 , NC) 
NBITE (6,699J) 

690J rORNAT (5X,*THE C NATRIX'/) 

IK (kS .LE. 10) GO TO 6970 
WRITE (6,6960) (HOG (L) ,L-1, JO) 

DU 6921 1-1, NT 

6921 WRITE (6,6901) kAN J (T) , (C (I , J) , J- 1, kCOL 1) 
WRITE (6,6960) (HDG(L) ,L-)1,60) 

DO 6922 1-1, NT 

6922 WRITE (6,6901) NAN J ( I ) , (C ( I , J) , J- NCOL2, NS) 
CO TO 6975 

6970 WRITE (6,6960) (HOG (L) , L- 1 , JO) 

DO 6929 1-1, NX 

6923 WHITE (6,6901) NAN 1 ( I ) , (C ( I , J) , J- 1 , NS) 

6975 WiUTb (6,6904) 


o n o n n n 


X) 


ORKilNAL HAUK 
OF FOUR gt’Ai.lTY 


690« P'»M»t(9X,'THE D 

k*m (6,6960| (AOaiLI ,L*1> iO) 

DO 6VJI t*1,NT 

69JI M9!Tr 16,6901) 669 1 ( l| , (0 (I , J» , J* 1 . «C) 

IP OPSIdC k PPR96BPKT DITASET WITH A.B,C,D ntrilCES 
I99LI-T 6RITL STPTS KCPE 6ITH OESIIBD PO>9«r 


••lErOF riASSAGE MUIblli PPIBTOUT 

IP (6660P .91. 0) «Pirr)6,711) HEEOI 
711 17U1kT (///,• ••• AT LEAST ‘.IJ,* PPPOIS 

STOP 
L7C 


PEOCDEP READIN 


nrnro'-p ►' 

PAT AM nn r 

IF ('"TTP'«"); ro TO 'FMP' 

ni *?"i AY 'F"Tn oata - Fii!i 511 (II Til _riir.' 

HA^r nriK 

FPi T m t' 
nFFAlMT SV.M'Y-'’ 

Trri T nnp 
ysFi I.'" nriP;o 
nFFAIll.T SVSI PX-r 

ni50|AV 'VT-PATA^FT rOMPLFTl.o - 1*001 FY IF *'EC' 
*E*in' 
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APPKNDIX C 
'QINFO' dataset 


Format : 


line 1 NS, NC. ^^Y 

NAME array 
SF array 
NAM2 array 
NAM3 ariay 
NDF.LT array 
NDRIV array 
NOBSV array 

last line NCOL array 


313 

30A4 

8(10F8.2) 

30A4 

30A4 

3013 

3013 

3013 

3013 


Descrlpt ion : 


NS : 

number of 

states 

NC : 

number of 

inputs 

NY : 

number of 

outputs 

NAME: 

names of 

the states and Inputs in the order perturbed for 


generation of hybrid data, 
length of each name - 4 alphanumeric characters 
Each position of the NAME array contains the name of a state 
or an input. Since each state and each input is perturbed 
in both a plus and minus direction, each position of the 
NAME array corresponds to two readings and hence, two dis- 
turbance numbers. The correspondence between a position In 
the NAME array and two disturbance numbers is as follows: 

N: position in NAME array 

(e.g., NAME(N) ■ STl name of state 1) 

disturbance number of plus perturbation 2*N 
disturbance number of minus perturbation 2*Ni-l 

NCOL: determines columnar structul'e of the matrices 

The integer values represent the positions of the NAME array 
in the order desired for the columns of matrices. 

(e.g., NAME array STl, ST2, INI, 1N2, ST3, ST4„ 1N3, IH4 
position 12345678 

Want columns 1, 2, 3, and 4 of A and C matrices to be 
with respect to STl, ST2, ST3, ST4 and want columns 1, 2, 

3, and 4 of B and D matrice to be with respect to 
INI, IN2, IN3, IN4. 

Then NCOL array is 1, 2, 5, 6. 3, 4, 7, 8) 


! 

I 

I 

4 


I 
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The data for the arrays SF, NDELT, NDRIV, und NOBSV depend on the order 

of the values In the rows of the DATA matrix. I.et 'ROW* mean the gen- 
eral form of a row In the DATA matrix. 

SF scale factors corresponding to values in 'ROW'. 

the real x values in 1/x scale factor representation. 

NDELT: positions in 'ROW' occupied by the values of the states and 

Inputs. 

NDRIV: positions in 'ROW' occupied by the values of the state deriv- 

atives ordered as the rows of A and B matrices are desired. 

NOBSV: positions in 'ROW' occupied by the values of tne outputs 

ordered as rows of C and D matrices are desired. 

NAM2: names of the Ttate derivatives in order specified in NDRIV. 

length per name - A alphanumeric characters. 

NAM3: names of the outputs In order specified in NOBSV. 

length per name - 4 alphanumeric characters. 


^ I 
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APFKNDIX D 

DATA PROCESSING PROCEDURE 

Note: This procedure is specialized for use with the NASA-Lewls 

Research Center EAI Pacer 100 computer system. 

TEKTRONIX 4'JlO TERMINAL INITIALIZATION 

1. turn terminal ON (power light will come on). 

2. turn 3h0 switch to the ON position. 

3. turn SEL switch to the OFF position. 

4. set BAUD switch for 1200. 

5. set NORM/TP-WR switch for NORM. 

6. set LOCAL/LINE switch for LINE. 

7. on the phone modem, moraentalrly depress the CONNECT button. 

(the CONNECT light will come on when line Is connected). 

8. hit 'BREAK' key on terminal. 

9. 'logon' to 380 as normally done. 

(check TERMINAL ID; it should be posted). 


TEKTRONIX 4922 FLOPPY DISK INITIALIZATION 

1. turn POWER switch ON. 

2. depress RESET button. 

3. turn PROMPT MODE switch OFF. 

4. select proper disk. 

a. check disk notebook. 

b. make sure disk moves freely In envelope. 

5. set DISK switch for disk being used. 

a. A If disk 1. 

b. B If disk 2. 

8. set TERMINAL switch (switch #1) for terminal being used. 

a. LEFT If left terminal. 

b. RIGHT If right terminal. 

7. set TRANSMISSION switch (switch 02 ) to the LEFT, 

(LEFT setting - data transmission stops at 'XOFF' character). 

8. set R-W switch to MANUAL READ. 
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DATA TRANSFER FROM FLOPPY DISK 

1. (JepreHii RESET button on dink. 

2. set beglnninR track address In TRACK ADDRESS. 

3. depress LA button on disk. 

4. type at terminal 

READIN 'DTRK' (CR). 
when the procdef starts, the message: 

'ENTER DATA - FINISH WITH - END* 
will print at the terminal. Data transfer will be automatic 
the user can determine this by the periodic "noise" the disk 
will make during this transfer. 

5. type at terminal : 

_END (CR) 

when the end of the data Is reached. 

This will end the EDIT dataset created by procdef READIN. 

It may he necessary to type _END (CR) more than once before 
the 360 accepts it. These unaccepted lines are entered into 
the dataset 'DTRK' and must be deleted through REDIT 
before using 'DTRK' as Input to the MATCALC program. 

'MATCALC program EXECUTION 

After creation of datasets 'DTRK' and "QINFO', the program 
can be executed by typing at the terminal: 

DDEF FT05F001.VS, 'DTRK' 

DDEF FTORFOOl.VS, 'QINFO' 

DDEF FT09F001,VS, 'QOUTPUT' 

MATCALC 

NOTES : 

1. When the message 'MATRIX DATA READ IN' etc. is printed 

at the terminal, depress RESET button on disk. 

2. Final printout at terminal will be the A, B, C, and D 

matrices 

3. If errors are detected, appropriate messages will be 

printed out, allowing the user opportunity *o make 
corrections. The procedure for checksum error cor- 
rection follows. 
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CHECKSirM ERROR PROCEniJRE 

If the read and calculated values of the checksum for a line of 
floppy disk data are not equal, a checksum error has occurred. 

An error message will print at the terminal, allowing the data 
line to be re-entered as followst 

1. set LOCAI./LINE switch to LOCAL on terminal. 

2. depress RESET button on disk. 

3. set beginning track address in TRACK ADDRESS. 

4. depress LA button. 

5. at the terminal, 

the 'CTRLQ' key will input one line of' floppy disk data 
to the terminal CRT screen. 

a. Lines of data should be Inputted until the proper 

disturbance reading data is reached. This is done 
by comparing the RUN number outpi'tted in the error 
message to the values in the line containing the 
disturbance number for each reading. 

b. Lines of data from that reading should then be 

Inputted, noting the checksum for each line. When 
the checksum equal to the PREVIOUS checksum 
(printed in the error message) is reached, the 
LOCAL/LINE switch should be set to LINE. 

c. The next data line should be Inputed by the 'CTRLQ' 

key. This line will be accepted by the 360 as the 
re-entered data line. If the checksum of this data 
line is on the next line, that line will automatically 
input to the 360. 
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TABLE I. - PERTURBATION CONTROL 
CARDS FOR 4-STATE, 4-INPUT, 
4-OUTPUT EXAMPLE 


45 

T5 

41 n;' 

ri 


T5 

45 

71 

41 


Mr 

a? 

M?’ 


TABLE II. - PERTURBATION ORDER OF STATE 
VARIABLES AND SYSTEM INPUTS FOR 4-STATE, 


4-INPUT, 4-OUTPUT EXAMPLE 


Parameter 

Console ADC Channel 

STl 

Primary 

1 

ST2 

Primary 

2 

INI 

Primary 

5 

IN2 

Primary 

6 

ST 3 

Secondary 

33 

ST4 

Secondary 

34 

IN3 

Secondary 

37 

IN4 

Secondary 

38 

Note: The 

perturbat Ion 

order of 


state variables and system In- 
puts per console Is arbitrary. 
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TABLE III. - READING NUMBERS FOR 4-STATE, 
♦-INPUT, 4-OUTPUT EXAMPLE 
RdR. # (■ disturbance #) 


1 

BASE RUN 

2 

STl + 

3 

STl - 

4 

ST2 > 

S 

ST2 - 

f> 

INI * 

7 

INI - 

8 

IN2 

9 

IN2 - 

10 

ST3 - 

11 

ST3 - 

12 

ST4 * 

13 

ST4 - 

14 

IN3 

13 

IN3 - 

lb 

IN4 -E 

17 

TN4 - 


Order of values per reading: 

STl, ST2, DRVl, DRV2, INI, IN2, OUTl, OLT2, ST3, ST4, DRV3, DRV4, 
IN3, IN4, OUT3, OUT4 

Note - STJ represents state variable Xj where J ■ 1 to 4 
INJ represents input variable Uj where ./ ■ 1 to 4 
OUTJ represents output variable Xj where J ■ 1 to 4 
DRVJ represents state derivative Xj wh^re J ■ 1 to 4 


) 


I 


1 


I 


I 


I 





■ 

H 

■ 

■ 

■ 

B 

■ 

■ 

■ 

w 

1 

1 


16 l«2 
24)79i 

24< UO 

0 

2 

2 


If M9 
24* 210 

24*9) 

•11) 

1 

1 


16 <74 
24SS9 1 

24)92 

104 

« 

4 


U 194 
24* 177 

2)011 

0 

> 

6 


1* 199 
24 ta7 

24 112 

0 

ii 

6 


14 He 

24114 

tOe 

24< )01 

) 

7 


U 19 7 
24)229 

24197 

- ij) 

H 

9 


16 iH4 
24.94) 

24)14 

0 

'* 

C 


14119 

24499) 

24il7 

0 

11 

10 


1l 111 

44 .172 

24)9* 

0 

1 1 

11 


1 * M9 

24 jl7 

0 

44 j 14)1 
1 1 

12 


1* H9 
24‘ lit 

24)17 

1 

1 1 

1 1 


U 1 ■» 
2<t --42 

24114 

0 

l4 

14 


If 199 

44)17 

0 

24i 94) 

0 


1 i 

1) 


16 He 
24«9)2 

24>17 

0 

H 

I 4 


It H9 
24 61 *4 

49)17 

0 

17 

17 


16 H9 

24317 

u 



Poor quaitty 



I 
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TULI IV. - fLPfn DIM data VDI 4-«TATI. A-tMTirt, ♦-OVT'^ n’AlTU 


0 t»H1 JiiiVt l«ia« 2V»AI i«AO« UAM 

C 1vl«} iWiT* MTJV i«ivt 2«»0I U«K 

0 1*MI jmli 1604* 24S«4 3A40J 144 1 

•40A UMl 14TV9 2&0I0 24A0 1 U«1 

4eJ tbial ^4)7* UJt) 74049 ik403 14«1 

0 1«*VA ;«)7S UI9A ^11)41 JUtOi 1441 

0 IbOMi ^447} 14 1*4 ^440} l44l 

k»J 1«1%1 'BU! UT*‘. «44«l J*^02 1441 

>B07 It lei JkODk H )A! 2«l))1 24401 tfvl 

U 1el4t ^k:?i ItlAe 24''«2 2*094 U41 

0 14 141 2417* 1414S 24441 24120 U41 

0 Ulil 24)7* 1)H) 24)41 24401 1474 

0 1cl4l .4)7) 1614) 24)91 24601 1610 

0 1oJ4J 24)7) 1614) 24)41 24401 1441 

0 16141 24)7) 14 Id) 24)60 2*401 1(4l 

0 16171 24)7) 1)16) 24)91 24601 14416 

0 14161 24)7* 1614) 24)61 24401 16414 


t 16 192 24)44 1-»Jt 

0 24)74 1414) 24)4) 1441 ' 

0 24)7t 16 It) 2«)9) lr* lu 

0 24)74 14 )4) 24)9) 1»4l<- 

0 249'/6 16)9) 24)9* Itil 

0 24)74 16)9) 24)9) l*:41 . 

0 24)76 16 191 24)9) 1t4ir 

0 24)74 16 U) 24)9) I^b1t 

0 24)76 16)4) 24)9) P41' 

0 24)76 14142 2)099 1441 

0 24)74 16)9] 24112 1441 

• 141 24)77 1f|91 24)9) It 7. I 

279 24)77 16)9) 24)96 1»09> 

0 2)049 14)61 24)4) I44lr. 

0 24096 1619) 24)9* iMilt 

276 24)77 U70f 24)9l 1t4lD 

■142 24)77 H0)9 24)9t 14411. 
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OHIorNAL FA<iR K' 
OF OUAI.ITV 


It 


TABLE V. - Musi' abbay and 'data* natrib rr>B •-btatb. 4>inttrr, A-ovmrr bia»«l( 

XBUi ribat 
1 
i 
i 
A 

s 

( 

7 

A 

» 

13 
II 
U 
I I 
1A 
1i 
It 
17 


DATA RArtll 

Coluan 


Rw 

1 

2 

3 A 

1 

6 

7 

8 

1 

1.».'10a9 

1.47t4) 

0.040003 3.000000 

0.999817 

1.87486 

1.00000 

1.874)4 

2 

1.L2043 

1. 37143 

•.204464 1*0 10. 000 000 

0.9999)9 

1.87492 

1.01484 

1.874)8 

1 

O.ttr 1201 

1.47622 

1. 14R*47E-011. 000000 

0.9999)9 

1.67492 

0.980469 

1.87461 

4 

1.< 3024 

1. )142« 

1.C00000 -.)873’»3!- 

010. 999939 

1.37492 

i.ooooe 

1.91 143 

> 

1.01024 

1. A I96 0 

0.000000 3. )5)241)t> 

010. 999919 

1.87492 

1.0000( 

1.6)461 

4 

1.i)0U24 

1.«7133 

0. 143347|(-010. 000000 

1.01917 

1.97492 

1.00006 

1.874)8 

7 

1.1;3131 

1 . oTloO 

•.20446lt>011. 000000 

0. 479638 

1.87492 

I.ooooe 

1.87318 

R 

1.00024 

1.47633 

0.000000 0.353241!- 

010.999939 

1.912)1 

1.00006 

1.873)8 

9 

1.1)0031 

1.47660 

0.000000 -. 3466101- 

010. 399939 

1.83746 

1.00O0( 

1.873 14 

10 

1.1)0043 

1. 37475 

0.000000 3.000000 

0.999617 

1.87492 

1.00000 

1.87346 

11 

1.00024 

1.47(60 

0.010000 0.000000 

0.999939 

1.87492 

1.00006 

1.67418 

12 

1.000)1 

1. 37660 

0.000000 0.000000 

0.49993) 

1.87492 

1.00006 

1.87518 

13 

1 .(‘0024 

1 . 07ej3 

1.010003 3.000000 

0.9999 )« 

1.47492 

1.00006 

1.C7314 

14 

1. 40024 

1.37(60 

0.000000 0.000000 

0. 9999)9 

1.87492 

1.00006 

1.874)8 

13 

1.01124 

1. m76o0 

0.000003 0.000000 

0. 9999 )9 

1.47492 

1.00006 

1.87511 

16 

I.LOOJI 

1.47(60 

0.000000 0.000000 

0.9999)9 

1.97492 

1.00006 

1.873)4 

17 

1.03024 

1.47660 

0.030009 1.000000 

0.9999)9 

1.8-492 

1.00006 

1.875)4 


9 

10 

11 12 

11 

16 

13 

16 

1 

3.7)436 

2.10427 

0.000003 0.000009 

3.75000 

1.99976 

). 73273 

2.00105 

2 

3. ■'34 12 

2.00)91 

0.000000 0.000000 

1.73000 

1.99988 

1.75290 

2.00)17 

3 

3. "Saw 

2.00)91 

O.OOnOO) 1.009000 

3.73000 

1.99984 

1. 75290 

2.00)17 

4 

3. :i4 12 

2.10)91 

9.000000 0.000000 

3.73000 

1.49988 

), 74290 

2.00 )17 

3 

3. •'3397 

2.00)91 

O.OOOCOO 0.000000 

1.73000 

1.94988 

1.74290 

2.00)17 

6 

3.73)97 

2.0C)')1 

0.000000 0.000000 

1.73000 

1.9948H 

3.732)0 

2.00)17 

7 

3. Via 12 

2.00391 

0.000009 9.000000 

1. 73000 

1.99988 

1. 73290 

2.09)17 

A 

3. •’3)97 

2.00311 

0.000000 0.000000 

1.73000 

1.99984 

1.73290 

2. '0317 

9 

3. '7412 

2. JU391 

0.000003 1.000000 

1. 73000 

1.99988 

). 74 2 9 0 

2.00)17 

10 

).' *965 

2.00413 

•.782776!>010. 000000 

1.73000 

1.99976 

).S2829 

2.00317 

11 

3.44042 

2.00)74 

0.700374e>019. 009000 

3. 73000 

1.49988 

1.67920 

a.r ,0317 

12 

3 .V -4 12 

2. J444 3 

3.000000 -.443113!- 

01 1.75013 

l.'998f 

1.73290 

2.04182 

13 

3.734 12 

1 . 3a543 

O.0C9003 3. )40376E- 

01 3.75013 

1.99J74 

1.79305 

1.96472 

14 

).Vt4 12 

2.00)11 

1.6)863)1-010.000000 

1.8232) 

1.99988 

3.73290 

2.00117 

13 

3. V3412 

2.00)31 

-.779724E-019. 000000 

). 67334 

1.99988 

1.73290 

2.00)17 

16 

3.7 412 

2.00)91 

0.000000 0.3)6914!- 

01 3.75013 

2.03931 

). 73)05 

2.00117 

1? 

3.V3412 

2.00)91 

0.000000 -.441893!- 

01 3.75015 

1.9603) 

3. 73)05 

2.00)17 


ORIGINAL PAGE Ih 
Cr POOR QUALmS 
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TALU VI. • 'qiNro' MTAsrr mt 4-«'i *i. 4-iarirr, 4-o>jrPvr exawu 

4 » 4 

STI St2 111 112 ST« 1*3 1*« 

2. 2.) 3. 3.) 3. 3.S 3. 3.S S. « 

«. «. S. 4. i. *. 

MV1t.NV3b*VJ0«'’4 
OOT1UUT3oar J0UT4 

1 3 ^ « 9 13 11 1 * 

J « 11 13 

7 4 IS 19 

1 / S 9 1 4 7 * 


Hocal Appandli C glvaa tha raqulrad foraat apaclf teat Iona for 'QlNfO' dacaaat 
and dafinaa tha arrava In tha NATC/.IX p.ofraa which ara Inltlalliad by 
thla data. 


I 


OHK^INAL b 

OF PO)R quality. 
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TAIU Vlt. - 

A. A » C f 

D HATRIX HESULTI POP 

4-4TATE , 




4 INPUT. 

4-oirmrr example 


THE 

A 

RATH : < 






S. I 

S?? 

in 

ST4 

D»V 1 


-0. V 1 H 

u.ooco 

0.000] 

O.UOOO 



O.ODU 

-O.KIrt 

3.0000 

0.0300 

DRV ( 


O.O'UU 

J.OOOO 

-0. V 01 

0.00.0 

DRV* 


u.ono 

0.0000 

0.0000 

•0. ‘■121 

Tiie 

J 

RATP ti 






rn 

IK2 

IVJ 

!N4 

OfcV 1 


0.1 

0.0000 

0.0000 

0 . cooo 

DRVi 


0.0 J 30 

u. vottO 

0. 0000 

O.UOOO 

OR / i 


o.oo.io 

n.nooo 

0. ».»7 d 

0.0000 

0 PV 4 


O.O'MO 

0.0000 

0. 0000 

0. >64 1 

THE 

J 

RATH :t 






i. 1 

5X1 

in 

STtt 

ou: 1 


1. ».)l 

u.oooo 

-0. 4 JIOP-OI 

O.COOO 

ou rv 


-C.5 M^t-12 

1. 002 

3.M12E-3) 

0. 0000 

ou:i 


O.O'UO 

J.OOCC 

O.'lOlO 

0. 11 J 2» 

our 4 


0.0 ) )0 

0.0000 

O.OOOO 

1.002 

THE 

0 

RATS :< 






1 1 1 

IN2 

IN I 

IR4 

OUT! 


0.0)0 u 

0.0000 

0..1000 

0. 0000 

ou;^ 


0.1) . V )0 

0.0000 

3. *0*75-01 

0.0000 

OUT! 


0.0 iin 

o.oono 

0. 3000 

0.0030 

OUiV 


).0) J 0 

o.ooou 

0.0300 

O.COOO 








